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1.0  INTRODUCTION 

Many  missiles  have  a  spherically  capped  blunt  nose.  For  missiles  designed  for  high 
angle  of  attack,  a  blunt  nose  is  sometimes  preferred  for  the  purpose  of  increasing  stability. 
Therefore,  a  hemisphere-cylinder  can  be  selected  as  a  first  step  toward  understanding  the 
flow  field  over  blunt-nose  bodies.  The  reasons  for  selecting  a  hemisphere-cylinder  are  as 
follows:  (1)  a  hemisphere-cylinder  is  simple  in  shape  and  can  be  used  in  the  approximation 
of  the  flow  field  at  the  nose  for  spherically  capped  bodies,  (2)  a  hemisphere  nose  retains 
its  identity  when  the  body  is  at  an  angle  of  attack,  and  (3)  the  hemisphere-cylinder  itself 
has  many  engineering  applications  such  as  the  well-known  use  of  a  flow  probe  or  sensing 
device  sh'ape.  In  this  report,  the  flow  field  about  a  hemisphere-cylinder  at  zero  incidence 
is  presented,  and  extension  to  an  angle  of  attack  will  be  reported  separately.  Understanding 
the  flow  field  about  a  hemisphere-cylinder  is  particularly  difficult  for  such  a  blunt  shape 
w'hen  the  free-stream  Mach  number  is  in  the  transonic  and  low  supersonic  flow  range. 
As  a  first  step  toward  understanding  the  flow  phenomena  about  a  hemisphere-cylinder, 
a  fundamental  and  systematic  study  of  the  flow  field  about  a  hemisphere  cylinder  at 
zero  incidence  is  presented  in  this  report.  Both  theoretical  and  experimental  results  are 
presented  for  the  Mach  number  range  from  0.7  to  2.0. 

As  shown  in  Fig.  1,  the  flow  past  a  hemisphere-cylinder  may  be  described  for  two 
different  cases,  M^  >  1  and  M^  <  I .  For  M^  >  1 ,  the  flow  field  influenced  by  introducing  the 


b.  M  <  1 

eo 

Figure  1.  Description  of  flow  field  about  hemisphere  cylinder. 
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hemisphere-cylinder  is  restricted  to  the  space  enclosed  by  the  bow  shock  and  the  body 
as  shown  in  Fig.  la.  Behind  the  bow  shock,  the  flow  is  divided  by  a  sonic  line  into 
a  subsonic  region  at  the  nose  portion  and  a  supersonic  region  downstream.  On  the  other 
hand,  for  <  1,  the  influence  of  the  body  covers  all  the  flow  space.  As  Mach  number 
increases,  two  flow  phenomena  appear:  a  supersonic  pocket  and  a  normal  shock  extending 
from  the  body  as  shown  in  Fig.  1  b.  As  -»■  1  subsonically.  the  supersonic  region  increases 
and  the  normal  shock  moves  downstream.  In  principle,  the  bow'  shock  moves  upstream 
to  infinity  if  =  1  is  approached  supersonically,  and  the  normal  shock  moves  to  infinity 
downstream  (provided  the  cylinder  is  infinitely  long)  as  M^=  1  is  approached  subsonically. 
Consequently,  there  is  no  shock  in  the  flow  field  at  M_  =  1 .  However,  in  reality,  the 
cylinder  is  finite,  viscous  effects  exist,  and  shock  systems  will  exist  at  =  1.  It  is  the 
purpose  of  this  report  to  study  how  well  these  flow  phenomena  can  be  predicted  by 
inviscid  theory  and  to  evaluate  the  accuracy  of  inviscid  theory  by  comparison  with 
experimental  data. 

Numerous  theoretical  studies  of  inviscid  supersonic  flow  past  blunt-nose  bodies  are 
reported  in  the  literature  (Refs.  1  through  6);  however,  the  Mach  numbers  considered 
are  generally  greater  than  2.  Only  a  few  studies  (Refs.  2,  3,  5,  and  6)  for  spherical  nose 
shapes  have  been  carried  to  Mach  number  1,5  or  1.2,  and  the  computations  were  limited 
to  just  a  portion  of  the  hemisphere.  Among  the  various  calculation  techniques  used,  many 
investigators  (Refs.  1  through  4)  prefer  the  direct  method  of  a  time-dependent, 
finite-difference  solution  to  the  unsteady  Euler's  equations.  The  steady  solution  is  obtained 
asymptotically  as  time  advances.  Even  though  this  scheme  has  shown  satisfactory  results 
for  a  spherical  nose  for  >  1.5.  the  applicability  of  the  method  has  not  been  investigated 
for  M^  <  1.5.  In  this  report,  the  numerical  methods  of  Refs.  2  (modified  and  referred 
to  as  Code  A)  and  4*  (referred  to  as  Code  B),  based  on  the  time-dependent  finite-difference 
solution,  were  examined  to  check  the  validity  of  the  solution  at  supersonic  Mach  numbers 
as  low  as  1 .048  for  the  hemispherical  nose  portion.  For  the  afterbody,  the  time-dependent 
scheme  as  described  in  Refs.  2  and  4  is  not  applicable  at  M^  <  2;  therefore,  the  method 
of  characteristics  (Refs.  7  and  8)  was  utilized  for  the  computation  since  the  flow  is 
supersonic  everywhere  downstream  of  the  tangent  points. 

Not  many  calculations  have  been  reported  for  blunt-nose  bodies  of  revolution  in 
subsonic,  supercritical  flow.  To  the  author's  knowledge,  only  two  papers  present  solutions 
to  this  type  of  problem.  In  Ref.  9,  a  time-asymptotic,  finite-difference  method  for  solving 
the  nonisentropic  Euler's  equations  and  exact  boundary  conditions  for  closed  blunt  bodies 
is  presented;  in  Ref.  10  a  finite-difference  relaxation  method  is  presented  for  numerical 
solution  of  the  full  potential  equation  and  exact  boundary  conditions  for  general 


*The  computer  program  was  made  available  to  AEDC  from  Mr.  Putz  of  Sandia  Laboratory. 
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axisymmetrical  bodies  (blunt  or  pointed  nose).  The  method  of  Ref.  10  is  particularly 
convenient  because  it  can  treat  bodies  witli  an  open  end  and  is  also  applicable  to  a 
supersonic  free  stream.  Therefore,  it  has  been  chosen  for  the  hemisphere-cylinder  study 
presented  in  this  paper  and  is  referred  to  as  Code  C.  It  is  well  know'ii  that  in  subsonic, 
supercritical  flow,  the  interaction  of  the  shock  and  boundary  layer  strongly  changes  the 
tlow  as  predicted  by  the  inviscid  theory.  The  effects  of  the  interaction  for  the 
hemispliere-cylinder  studied  w'ill  be  discussed  later  in  the  presentation  of  results. 

Experiments  of  supersonic  flow  past  a  spherical  nose  as  related  to  the  present  study 
have  been  reported  (Refs.  11  through  13).  Data  are  available  for  the  shock  standoff  distance 
and  bow  shock  shape  for  =  1.048  to  1.82.  However,  surface  pressure  data  for  a 
hemisphere-cylinder  are  available  for  >  1.42  only.  Shadowgraph  and  pressure  data 
to  evaluate  the  theoretical  calculations  were  obtained- at  Mach  numbers  from  0.7  to  1.3 
in  the  AEDC  Propulsion  Wind  Tunnel  (IT)  for  a  l-in.-diam.  lO-in.-Iong  hemisphere-cylinder. 

2.0  THEORETICAL  STUDY 

Three  different  computational  schemes  have  been  employed  in  the  present  theoretical 
study  of  an  inviscid,  compressible  tlow  past  a  hemisphere-cylinder  at  zero  incidence  for 
the  Mach  number  range  from  0.7  to  2.0.  The  reason  is  simply  that  there  is  not  a  single 
method  that  can  be  effectively  applied  to  cover  the  entire  flow  range  to  be  calculated. 

2.1  TIME-DEPENDENT  FINITE-DIFFERENCE  SOLUTION  TO 

THE  UNSTEADY  EULER'S  EQUATIONS 

This  method  is  applied  to  the  hemispherical  nose  portion  as  shown  in  Fig.  2a  for 
M^  >  1.  It  is  assumed  that  the  sonic  line  is  located  ahead  of  the  tangent  point  (this 

assumption  holds  for  M^  >  1.05.  approximately:  see  Appendix  A),  and  so  is  the  line 

of  limiting  characteristics.  As  a  result,  the  flow  downstream  of  the  tangent  point  is 
supersonic  everywhere  and  has  no  influence  on  the  tlow  upstream  of  the  tangent  point. 
It  is  well  known  that  the  partial  differential  equations  for  the  steady  How  in  this  region 
are  of  mixed  type  (i.e.,  elliptic  in  the  subsonic  region  and  hyperbolic  in  the  supersonic 
region)  and  are  difficult  to  solve.  One  way  to  cope  with  the  mixed  elliptic-hyperbolic 
type  of  differential  equations  is  to  use  the  unsteady  flow  equations,  which  are  always 
of  the  hyperbolic  type.  The  steady  solution  is  then  obtained  by  advancing  the  time. 

Among  the  many  techniques  reported  for  applying  the  time-dependent  method  to 

blunt-nose  bodies,  the  direct  method  (Refs.  1  through  4)  was  chosen  for  this  study.  The 

computer  programs  (Refs.  2  and  4)  with  necessary  modifications  were  used  in  the  present 
calculation  of  the  flow  field  about  the  nose  portion.  A  brief  description  of  the  equations 
and  numerical  scheme  used  in  obtaining  solutions  is  given  in  the  following  paragraphs. 
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a.  Time-dependent  finite-difference  metiiod 


r 


COMPUTATION  GRID 


LAST  STATION 
USED  IN 
COMPUTATION 


c.  Relaxation  method 
Figure  2.  Methods  of  calculation. 
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2.1.1  Differential  Equations 

The  system  of  unsteady  Euler  equations  written  in  a  body-oriented  coordinate  system 
(s,n)  arc  as  follows; 


_  (9p  d(rpa) 

^  (9t  ^  ds 

Krpv  =  0 

(1) 

dll  du 

K  —  +  11  —  + 

dl  ds 

du  1  (9p  „ 

KV  -r-  +  -  =  K.UV 

dn  p  ds 

(2) 

dv  dv 

U  -IT-  + 
dt  dS 

dv  K  dP  ^2 

KV  =  -K  U 

dn  p  dn 

(3) 

5s  •  ds 

K  3-  +  U  3-  - 

01  ds 

ds 

KV  =  0 

dn 

(4) 

P  =  pye-^'S- 

(5) 

It  should  be  noted  that  Eq.  (1)  is  singular  at  s  =  0  and  hence  is  replaced  by 


K 


2^ 

dl 


„  dip  ii)  diK^pv) 
ds  ^  <9n 


(6) 


In  Ref.  4,  both  spherical  and  cylindrical  coordinate  systems  were  used.  In  general,  the 
system  of  equations  is  nondimensionalized  and  transformed  into  the  rectangular  calculation 
grid  as  shown  in  Fig.  2a. 


2.1.2  Numerical  Schemes 

1.  For  the  interior  points  in  the  computation  grid,  Ref.  2  approximates  the 
spatial  partial  derivatives  by  a  central  difference  and  the  time-partial 
derivatives  by  a  forward  difference.  Stabilizing  terms  for  artificial  viscosity 
are  added  in  such  a  way  that  the  effects  of  the  terms  were  made  negligibly 
small  (see  Section  4.0).  In  Ref.  4,  the  interior  points  are  computed  by 
the  implicit  predictor-corrector  scheme  of  Ref.  14. 

2.  The  bow  shock  boundary  conditions  in  both  schemes  are  represented  as 
a  discontinuity,  with  the  iteration  scheme  suggested  in  Ref.  1 . 

3.  The  body  boundary  conditions  in  Ref.  2  are  represented  by  backward 
difference  in  the  normal  coordinates  where  v  =  0  on  the  body  surface  (the 
other  differences  remain  the  same  as  for  the  interior  points),  whereas  in 
Ref.  4  the  body  streamline  conditions  are  satisfied  by  the  same  iteration 
scheme  used  in  the  shock  boundary  conditions. 
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In  both  Codes  A  and  B  the  time  interval  that  satisfies  the  stability  requirement  is 
automatically  built  in;  therefore,  one  needs  only  to  specify  the  nodal  points  and  the  number 
of  advances  in  time.  The  steady  condition  is  reached  when  the  computed  bow  shock 
velocity,  Ws,  on  the  axis  is  sufficiently  small  (on  the  order  of  lO"'*  to  lO'^  of  the  free-stream 
sound  speed). 

2.2  METHOD  OF  CHARACTERISTICS 

The  flow  field  about  the  cylinder  portion  was  computed  by  the  method  of 
characteristics.  The  General  Applied  Science  Laboratories  (GASL)  three-dimensional 
method-of-characteristics  code  (Ref.  7)  and  the  Lockheed  axial  symmetric  internal  flow 
code  (Ref.  8)  were  used  in  the  present  study. 

The  characteristic  equations  in  both  codes  were  derived  in  detail  in  Refs.  8  and  15 
and  hence  are  not  given  here.  The  calculation  proceeds  from  the  initial  station  toward 
the  downstream  stations  at  each  point  as  shown  in  Fig.  2b.  It  should  be  mentioned  that 
in  applying  the  method  of  Ref.  8,  which  was  developed  for  nozzle  flow,  a  cylindrical 
wall  is  artificially  put  in  with  a  radius  equal  to  the  bow  shock  distance  at  the  initial 
station,  i.e.,  the  tangent  point  (see  Fig.  2b).  As  a  result,  as  increases,  the  bow  shock 
moves  closer  and  closer  to  the  body,  and  the  application  of  the  Lockheed  code  deteriorates 
as  calculation  proceeds  downstream.  Therefore,  in  the  present  study,  the  Lockheed  code 
was  used  only  for  <  1.3  (where  the  artificial  wall  is  sufficiently  far  from  the  body 
and  has  insignificant  effects).  The  results  obtained  for  the  same  initial  data  provide  a 
comparison  of  the  different  computational  schemes  used  by  these  two  codes. 

2.3  RELAXATION  SOLUTION  TO  THE  STEADY  FULL  POTENTIAL  EQUATION 

For  <  1 ,  the  partial  differential  equations  are  elliptic  in  nature;  therefore,  the 

entire  flow  field  from  the  hemispherical  cylinder  to  infinite  space  must  be  computed. 
When  <  1,  the  assumption  of  irrotational  flow  is  approximately  valid.  Hence  a  solution 
to  the  steady  full  potential  equation  will  describe  the  flow  field.  It  is  given  as  follows 
(Ref.  10): 
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The  relation  between  ^  and  the  velocity  components  is 


u 


cos  6  + 


1  ^ 

K  ds 


V 


=  -sin  6 


diU 


Along  the  axis,  Eq.  (7)  reduces  to  its  limiting  form. 


2 

k“ 


+ 


0 


The  boundary  conditions  to  be  satisfied  are 

i/f  -•  0  as  n  -*  oo 


(8a) 

(8b) 


(9) 


(10a) 


and 


0  or 

an 


sin  0  on  the  body  surface 


(10b) 


As  described  in  Ref.  10,  sheared  cylindrical  coordinates  are  used  for  the  afterbody.  For 
the  hemisphere-cylinder,  Eqs.  (7)  through  (10b)  remain  the  same  for  the  sheared  cylindrical 
coordinates.  As  shown  in  Fig.  2c,  by  coordinate  transformation  the  flow  field  was  mapped 
into  a  rectangular  computation  plane  (^n),  where  the  relaxation  schemes  were  carried 
out.  For  >  1,  infinity  is, located  about  one  body  length  behond  the  body  end  station, 
whereas  for  >  1,  the  end  station  is  used  because  the  flow  is  supersonic  there. 

3.0  EXPERIMENTAL  INVESTIGATION 

Experiments  were  carried  out  for  the  purpose  of  evaluating  the  theoretical 
computations.  As  shown  in  Fig.  3.  a  hemisphere-cylinder  1  in.  in  diameter  and  10  in. 
long  was  sting  mounted  in  the  AEDC  Transonic  Wind  Tunnel  (IT).  A  description  of  the 
wind  tunnel  may  be  found  in  Ref.  1 6.  Pressure  data  were  taken  for  Mach  numbers  from 
0.7  to  1.3.  Shadowgraphs  were  taken  for  M^  =  0.8,  0.9,  0.95,  1.0,  1.05,  1.1,  1.2,  and 
1.3.  The  shock  standoff  distance  and  shock  shape  for  M_  >  1  and  the  normal  shock 
on  the  cylinder  for  <  1  were  obtained  from  the  shadowgraphs. 

In  addition  to  the  experimental  data  described  above,  the  experimental  results  of 
Refs.  1 1  through  1 3  were  also  used  to  evaluate  the  theoretical  results. 
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Figure  3.  Model  and  testing  setup. 


4.0  DISCUSSION  OF  RESULTS 

In  the  following  paragraphs,  the  theoretical  results  are  compared  with  other  theoretical 
calculations  reported,  in  the  literature  and  with  experimental  data  available.  The  validity 
of  each  method  of  calculation  is  discussed  and  evaluated.  To  facilitate  presentation,  results 
are  described  in  two  parts,  >  1  and  <  1. 

4.1  LOW  SUPERSONIC  MACH  NUMBER  RANGE,  2  >  M.  >  1 

4.1.1  Accuracy  of  Computation  and  Computer  Time 

For  >  1,  the  flow  about  the  hemispherical  nose  portion  was  calculated  mainly 
by  Code  A.  Computation  with  Code  B  was  also  made  for  comparison.  Since  stabilizing 
terms  are  added  in  Code  A  and  Code  A  results  were  used  as  initial  values  for  the  continuing 
computation  of  the  cylinder  portion  by  the  method  of  characteristics,  the  accuracy  of 
Code  A  must  be  evaluated  first.  The  accuracy  is  judged  by  two  quantities,  the  stagnation 
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pressure  and  the  total  enthalpy.  For  all  cases  computed,  the  stagnation  pressure  was  off 
by  no  more  than  l.S  percent  and  the  deviation  of  total  enthalpy  was  no  more  than  1 
percent;  this  is  considered  to  be  satisfactory.  The  values  of  pressure,  P,  and  the  u  and 
V  components  at  the  end  station  of  the  hemispherical  nose  are  tabulated  in  Appendix 
A.  These  values  were  used  as  the  intial  conditions  in  the  computer  program  for  the  method 
of  characteristics. 

For  the  IBM  370-155  computer,  with  which  all  the  computations  were  performed, 
the  computation  time  required  about  3  x  lO'^  sec  per  nodal  point  per  time  step.  As 
Mach  number  decreases,  the  nodal  points  and  the  time  steps  required  for  good  results 
increase  rapWly.  The  approach  to  steady  state  is  judged  by  the  ratio  of  the  shock  speed, 
Wg,  to  the  free-stream  sound  speed.  A  list  of  the  mesh  size,  time  steps,  and  shock  speed 
for  some  typical  Mach  number  values  is  given  in  Table  1. 


Table  1.  Hemisphere  Nose  Computational  Data 


M 

as 

Nodal  Points 

Time  Steps 

Shock  Speed,  Ws/v/p^P„ 

At  Axis 

At  Shoulder 

1.82 

12  X  12 

1,610 

0.1  X  10-4 

0.7  X  10-3 

1.60 

12  X  12 

2,531 

0.1  X  10-4 

0.8  X  10-3 

1.50 

13  X  13 

3,500 

0.7  X  10-5 

0.6  X  10-3 

1.42 

16  X  16 

4,451 

0.1  X  l(y4 

0.9  X  10-3 

1.30 

16  X  16 

8,000 

0.2  X  10-4 

0.1  X  10-2 

1.20 

24  X  24 

12,000 

0.7  X  l(y4 

0.6  X  10-2 

1.10  ■ 

30  X  30 

22,000 

0.4  X  10-3 

0.3  X  10-1 

1.05 

30  X  30 

30,000* 

0.14  X  10-2 

0.6  X  10-1 

'Discontinued  because  of  too-dow  convergence. 


For  example,  at  M^  =  1.5,  a  mesh  of  13  x  13  and  3500  time  steps  gives  good  results; 
however,  for  M^  =  1.2,  one  needs  24  x  24  mesh  and  12,000  time  steps.  It  is  important 
to  point  out  that  the  surface  pressure  approaches  steady  state  in  much  fewer  time  steps 
than  does  the  shock  standoff  distance.  For  example,  at  M____  =  I.l  with  2000  time  steps 
the  pressure  reaches  its  final  value  within  1  percent  at  the  nosetip  and  6  percent  at  the 
tangent  point,  and  with  4000  time  steps  the  corresponding  values  reduce  to  0.5  and  4 
percent,  respectively.  In  contrast,  the  shock  standoff  distance  at  the  nose  is  34  and  14 
percent  of  the  steady-state  value  with  time  steps  of  2000  and  4000,  respectively.  Figure 
4  shows  the  asymptotical  approach  to  the  steady-state  shock  standoff  distance  for  M^ 
=  1.2,  1.1,  and  1.05.  It  should  be  noted  that  for  M^  =  1.05,  the  computation  stops  at 
30,000  time  steps  because  of  too  slow  a  convergence,  and  the  steady-state  A/R  value 
was  obtained  by  extrapolation. 
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Figure  4.  Shock  standoff  distance  versus  time  steps. 

4.1.2  Shock  Standoff  Distance 

In  Fig.  5,  the  bow  shock  standoff  distance  A(R  is  plotted  as  a  function  of  free-stream 
Mach  number,  >  1.  The  curves  and  the  solid  symbols  are  from  various  theoretical 
results,  and  the  open  symbols  are  experimental  data  from  various  sources.  It  should  be 
noted  that  the  experimental  data  (Refs.  1 1  and  13)  are  for  spheres.  Since  the  flow  becomes 
supersonic  before  it  reaches  the  tangent  point  for  all  data  presented  in  Fig.  3  (see  Appendix 
A),  the  standoff  distance  and  the  bow  shock  shape  for  the  hemispherical  nose  should 
be  valid  for  comparison  with  the  sphere.  (Note  that  the  influence  of  limiting  characteristics 
at  1.0  <  <  1.1  remains  to  be  investigated.)  Two  values  of  7  =  1.4  and  1.13  are 

shown  in  the  figure  to  compare  with  the  available  data.  The  agreement  for  different 
theoretical  predictions  is  excellent  except  for  Ref.  6,  which  predicts  a  smaller  standoff 
distance  as  the  Mach  number  decreases,  and  Code  C  (Ref.  10),  which  consistently  predicts 
a  larger  value.  The  comparison  of  the  present  theoretical  results  with  experiments  is  also 
satisfactory  except  in  the  very  low  Mach  number  range,  1.05  <  M^  <  1.13,  where  the 
experimental  data  (Ref.  II)  are  lower  than  the  theory.  The  AEDC  experimental  data  for 
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air,  7  =  1.4,  are  consistently  less.  This  remains  to  be  investigated.*  Generally  speaking,  the 
shock  standoff  distance  for  a  hemisphere-cylinder  can  be  predicted  by  inviscid  theory  at  low 
supersonic  Mach  numbers. 


EXPERIMENT  THEORY 


□ 

AEOC,  PRESENT 

X-  1.4  ^ 

NASA  TR-RI,  NCF.  S 

CODE  A.  y  *  1. 4  1 

CODE  A.y  .I.Isl'^F  2 V RESENT 

r*  I.IS  f’STILP,  REF.  II 

A 

O 

y-  i.is  J 

k 

eiLINSKII  AND  LEBEDEV,  REF.  S 

0 

N«C«  TN-2000,  REF.  IS 

CODE  C  IREF.  10),  PRESENT 

or 

NPL.  REF.  12 

▼ 

CODE  B  IREF.  41,  PRESENT 

Figure  5.  Comparison  of  shock  standoff  distance  between 
calculation  and  experiment 


*See  Appendix  B. 
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4.1.3  Shock  Characteristics 


Different  methods  of  predicting  the  bow  shock  are  shown  in  Fig.  6.  The  agreement 
is  satisfactory.  Reference  6  gives  a  smaller  shock  layer  than  the  other  methods  as  Mach 
number  decreases.  Code  A  predicts  a  shock  shape  with  less  curvature  than  Code  B  for 
all  cases  studied.  Figure  7  shows  a  comparison  between  the  theoretical  calculation  and 
available  experimental  data.  The  agreement  between  the  present  calculations  and  the  data 


Figure  6.  Comparison  of  theoretical  shock  position. 
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of  Ref.  9  is  surprisingly  good  to  a  Mach  number  as  low  as  1.17.  Again  the  AEDC  data 
show  a  thinner  shock  layer  (see  Fig.  B-3  for  new  data). 


THEORY  -  CODE  A  (REF  PRESENT 


EXPERIMENT 


Mm 

M. 

0 

1.171 

1 

• 

1.05 
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1.30 

NACA  TN-2000. 
REF.  13 

4 

1. 10 
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1.62 

■ 

1.20 

1.81  . 

A 
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NPL, 

REF  12 
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Figure  7.  Comparison  of  theoretical  shock  position  witfi  experiment,  7  =  1.4. 
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4.1.4  Sonic  Lines 


The  sonic  lines  calculated  by  different  methods  are  compared  in  Fig.  8.  Calculations 
from  Codes  A  and  B  compare  well  with  each  other  and  with  results  of  Refs.  3  and  5 
for  >  1.5.  For  <  1.3,  the  sonic  line  predicted  by  each  method  starts  to  deviate 

M. 

'  All}  ^  '’•’ESENT 

O  I.S1 

O  1.3  GILINSKII  AND  LEBEDEV,  REE  5 
THEORY  Q  I.ZJ 

■  I.S  MORRETTI,  REF.  3 

-  CODE  A  (REF.  2),  PRESENT 

, - CODE  B  (REF  4),  PRESENT 


Figure  8.  Comparison  of  sonic  lines. 
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as  the  distance  away  from  the  body  increases.  The  relaxation  method  predicts  a  too-large 
subsonic  region  with  the  bow  shock  front  too  far  away  from  the  nose.  However,  in  the 
region  close  to  the  body,  the  sonic  lines  predicted  by  each  method  converge.  At  = 
1.2,  Code  B  gives  irregular  numbers  for  ail  the  flow  variables  in  the  re^on  close  to  the 
body  (see  Fig.  9  also).  Nevertheless,  the  shock  shapes  are  still  good,  and  the  part  of  the 
sonic  line  close  to  the  shock  can  still  be  sketched  as  shown  in  Fig.  8.  Significant  differences 
in  the  sonic  line  -location  are  shown  between  the  present  calculation  and  Ref.  5  for 
=  1.2.  At  the  present  time,  no  experimental  data  are  available  for  comparison  of  the 
sonic  line  location  and  hence  no  evaluation  can  be  made  of  the  prediction  methods. 


Figure  9.  Comparison  of  dieoretical  pressure  over  the  hemisphere  nose. 


4.1.5  Surface  Pressure 

Figure  9  shows  the  calculated  pressure  distribution  on  the  hemisphere  part  from 
various  sources.  The  agreement  is  quite  satisfactory.  As  described  previously,  at  =  1 .2 
the  results  from  Code  B  start  to  oscillate  about  Z/R  >  0.4  (the  same  mesh  number  was 
used  as  the  Code  A,  and  it  was  also  found  that  increasing  the  nodal  points  would  not 
help  to  improve  the  result),  and  it  will  not  run  for  <  1.2.  Thus,  for  <  1.2,  Code 
A  is  the  only  satisfactory  code. 

The  favorable  comparison  of  Fig.  9  for  the  nose  portion  gives  some  confidence  to 
use  the  data  obtained  from  Code  A  at  the  end  station  of  the  hemispherical  nose  to  start 
the  computation  of  the  flow  about  the  cylindrical  portion  by  the  method  of  characteristics. 
The  initial  data  of  P,u,v  used  in  the  present  calculations  are  shown  in  Fig.  10  for  (Y-l)/R 
<3.8.  A  test  run  was  carried  out  for  M^=  1.2  with  both  the  GASL  code  and  the  Lockheed 
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Figure  10.  Initial  data  for  MOC  calculation. 


code  with  the  same  initial  data.  It  was  found  that  the  GASL  code  requires  twice  the 
initial  nodal  points  (initial  data  were  linearly  interpolated  for  the  additional  nodal  points) 
as  the  Lockheed  code  to  obtain  the  same  accuracy  as  shown  in  Fig.  1 1 .  On  the  other 
hand,  the  Lockheed  code  is  unsatisfactory  for  the  downstream  surface  pressure  when 
>  1.3  (see  Section  2.0  for  explanation).  For  the  low  supersonic  Mach  numbers,  the  GASL 
code  was  used  for  >  1.1  (it  would  not  run  for  =  1.05).  and  the  Lockheed  code 
was  used  for  =  1.05  to  1.3.  The  overlapping  Mach  number  range  (M^  =  1.1,  1.2, 
and  1 .3)  provided  a  comparison  of  results.  Figure  1 1  shows  that  the  comparison  is 
satisfactory.  In  Fig.  12,  the  calculated  surface  pressures  are  compared  with  experimental 
data.  It  is  seen  that  for  all  cases  calculated,  the  time-dependent  method  plus  the  method 
of  characteristics  gives  very  good  agreement  with  experiment.  It  should  be  noted  that 
for  =  1.05  the  calculation  broke  down  in  the  cylinder  portion  at  Z/R  =  2.5  because 
of  too-low  Mach  number.  Calculations  were  also  carried  out  for  Mach  numbers  up  to 
1.42  by  using  the  relaxation  method  of  Ref.  10  (with  IBM  370-155,  a  complete  case 
using  a  mesh  halving  scheme  of  25  x  25,  49  x  49,  and  97  x  97  requires  approximately 
35  minutes  of  computer  time).  Surprisingly  enough,  the  agreement  between  the  relaxation 
method  and  the  experiments  is  good  up  to  M_^  =  1 .3,  particularly  for  the  minimum  pressure 
at  the  tangent  point.  For  M^  =  1.42,  the  pressure  predicted  by  the  relaxation  method 
is  not  satisfactory  and  hence  is  not  shown. 
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Figure  11.  Comparison  of  MOC  calculations. 


4.1.6  Oblique  Shock 

The  shadowgraphs  (Fig.  13)  show  that  systems  of  weak  oblique  shocks  exist 
downstream  of  the  tangent  point;  these  indicate  an  overexpansion  of  the  flow  for  most 
Mach  numbers.  The  major  or  first  oblique  shocks  each  have  nearly  the  same  shock  angle 
(i.e.,  about  40  to  50  deg  for  the  range  of  Mach  number  shown).  It  is  interesting  to  note 
that  in  plotting  the  characteristic  lines  obtained  from  the  Lockheed  code,  a  coalescing 
of  the  characteristic  lines  indeed  indicates  the  appearance  of  an  oblique  shock  as  shown 
in  Fig.  14  for  M^  =  1.2.  If  the  results  of  Fig.  13  are  compared  with  those  of  Fig.  14 
for  M^  =  1.2,  the  slope  and  location  of  the  major  oblique  shock  are  roughly  predicted 
by  in  viscid  theory. 

4.1.7  Flow  Characteristics  in  the  Shock  Layer 

To  describe  the  entire  flow  within  the  shock  layer,  a  mapping  of  the  pressure  and 
Mach  number  distribution  obtained  by  Code  A  and  GASL  code  for  M^  =  1.1,  1.122, 
1.2,  1.3,  and  1.42  is  shown  in  Figs.  15  and  16,  respectively. 
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Figure  13.  Shadowgraphs  of  flow  past  hemisphere-cylinder. 
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Figure  16.  Continued. 
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e.  M„  =  1.42 
Figure  16.  Concluded. 


4.2  TRANSONIC  MACH  NUMBER  RANGE,  0.7  <  M  <  1 

4.2.1  Supersonic  Pocket  and  Normal  Shock 

In  Fig.  17,  a  comparison  is  presented  of  the  location  of  the  supersonic  pocket  and 
shock  calculated  by  the  relaxation  method  of  Ref.  10  (note  that  the  method  of  Ref. 
10  does  not  use  the  conservative  differencing  with  a  "shock  point"  operator  as  does  Ref. 
17)  and  obtained  from  the  shadowgraphs  of  Fig.  13.  It  is  seen  that  the  shape  of  the 
shock  for  the  cases  shown  is  predicted  correctly,  but  the  location  is  shifted,  the 
experimental  results  show  the  shock  occurred  some  distance  downstream. 

4.2.2  Surface  Pressure 

In  Fig.  18  the  surface  pressure  calculated  by  the  relaxation  method  is  compared  with 
the  experimental  data  for  <  1 .  It  is  seen  that  the  theory  satisfactorily  predicts  the 
pressure  distribution  up  to  the  location  where  the  shock  seems  to  appear.  For  =  0.8 
there  is  stronger  disagreement  in  the  theoretical  and  experimental  pressure  distribution. 
The  corresponding  shadowgraph  given  in  Fig.  10  does  indicate  a  pronounced  boundary-layer 
separatation  for  =  0.8.  The  same  phenomenon  occurs  to  a  lesser  degree  at  other  Mach 
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Figure  17.  Comparison  of  supersonic  pocket  and  normal  shock. 


numbers.  Therefore,  it  is  assumed  that  the  difference  in  experimental  and  theoretical 
pressure  distribution  for  0.7  <  <  1.0  is  caused  by  the  interaction  of  the  shock  and 

the  boundary  layer. 

4.2.3  Sonic  Flow,  =  1 

Sonic  flow  presents  a  fascinating  case  in  both  theory  and  experiment.  One  is  pleased 
to  see  the  good  agreement  between  prediction  and  data  in  the  surface  pressure.  It  would 
be  interesting  to  point  out  that  the  disturbances  (not  shocks)  upstream  of  the  cylinder 
as  shown  in  Fig.  13  for  =  1  have  the  counterpart  in  the  numerical  computation  (Fig. 
19).  The  theory  indeed  predicts  a  supersonic  pocket  which  extends  to  infinity.  Also,  the 
sonic  line  moves  beyond  the  plane  of  the  shoulder.  With  finite  length  of  the  cylinder, 
the  shock  is  located  near  the  end  of  the  cylinder  as  predicted  by  the  theory.  Unfortunately, 
the  shadowgraph  for  =  1  does  not  cover  the  tail  end  of  the  cylinder  to  offer  a 
comparison. 
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Figure  18.  Comparison  of  calculated  pressure  with  experiment. 
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Figure  19.  Calculated  sonic  line  for  =  1  by  method  of  Ref.  10. 

5.0  CONCLUSIONS 

From  the  above  study  of  the  flow  field  about  a  hemisphere-cylinder  at  zero  incidence 
in  the  transonic  and  low  supersonic  Mach  number  range,  the  following  conclusions  can 
be  drawn: 

1.  The  inviscid  theory  can  satisfactorily  predict  the  flow  field  for  2  > 

>  l.OS.  The  method  used  is  a  combination  of  time-dependent  solution  for 
the  hemispherical  nose  portion  and  the  method  of  characteristics  for  the 
cylindrical  portion. 

2.  The  relaxation  method  of  Ref.  10  can  satisfactorily  predict  the  surface 
pressure  for  Mach  numbers  from  0.95  to  1.3,  but  the  external  flow  field 
predicted  is  not  as  good  as  the  surface  pressure.  For  Mach  numbers  near 
0.8,  experimental  results  show  an  interaction  of  the  shock  and  boundary 
layer;  hence  the  inviscid  theory  fails  to  give  satisfactory  results. 
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APPENDIX  A 

TABULATION  OF  FLOW  QUANTITIES  IN  THE  INITIAL  PLAN 
AT  TANGENT  POINT  ACROSS  THE  SHOCK  LAYER 


35 


M  =  1.048 

oo 

R 

U 

V 

RHO 

P 

M 

T 

i«OOOOOE  00 

1. 715006 

00 

070 

5.441006-01 

"4.350006-01 

1.621046 

■00“ 

/•99%CI6C:^TJ1~ 

1.171^1E  00 

1.567646 

00 

1.074936-01 

6.892066-01 

5.975006-01 

1.426286 

00 

8.669406-01 

l.J424[E  Od 

1.400006  00 

1.919966^1 

133626-01 

7.60000E-^1~ 

1.235516 

00 

9.34394EH1 

1.51362E  00 

1.348526 

00 

1.846266-01 

8.743436-01 

8.350006-01 

1.177136 

00 

9.550026-01 

1.60ABJE:  Off 

1.29S00E 

00 

1.770006-01 

9.327496-01 

9.100006-01  “ 

1.118376 

00 

9.756116=06 

1.836036  00 

1.267856 

00 

1.663996-01 

9.558156-01 

9.425006-01 

1.088336 

00 

9.860706-01 

2.0272VE  00 

1.240006 

iro~ 

1.  560046  —01 

9. /83976-01 

■  9.750006-01 

“T. 058096 

"00“ 

9.965286-01 

2.19845E  00 

1.230096 

00 

1.  454406-01 

9.920256-01 

9.924996-01 

1.04661E 

00 

1.000486  00 

2.36965E  00 

220006 

00 

1.3  50056^03 - 

1.005556  00 

17010006  00 

1.035096 

00 

: 17004436  00 

2.54086E  00 

1.213076 

00 

1.2644  76-01 

1.01667E  00 

1.025006  00 

1.024086 

00 

1.008196  00 

2.71206E  00 

1.200006  00 

1.  17996E-01 

17027716  00 

1.04000600 

1.013046 

00 

1*011966  00 

2.883276  00 

1.19652 6 

00 

1.114806-01 

1.031226  00 

1.045006  00 

1.008916 

00 

1.013366  00 

3.03448E  00 

1.193006 

00 

1.069966-01 

1.034726  00 

1.050006  00 - 

1.004786 

■00“ 

- 1.014766  00 

3.225686  00 

1.  192916 

00 

9.9995  76-02 

1.037006  00 

1.052506  00 

1.004266 

00 

1.014946  00 

1.396896  00 

1.192806 

00 

■9.499476^02 

1.039286  00 

17055006  00 

1.003736 

oir 

1.015126  00 

3.568106  00 

1.192806 

00 

9.24958E-02 

1.039516  00 

1.055306  00 

1.003546 

00 

1.015186  00 

T.  739306  TJO 

1.192806 

00 

8.999666^02 

1.039746  00 

1.055606  00 

1.003346 

00 

1.015256  00 

3.910516  00 

1.  192806 

00 

8.600126-02 

1.039646  00 

1.055606  00 

1.003056 

00 

1.015356  00 

4.081726 

"TJO — 

1. 192806 

T70 — 

1.039546 

■00“ 

- 1.055606 

■00“ 

- 1.002746 

“00 — 

““17015456 

DU 

4.252926 

00 

1.192816 

00 

7.  8003  56-  02 

1.039506 

00 

1.055656 

00 

1.002506 

00 

1.015546 

00 

■47424136 

00 

1.I92806“ 

■w 

7.400136-02 

1. 039466 

00 

1.0S570E 

00 

1.00224E 

00 

1.015436  DO 

4.595346 

00 

1. 192806 

00 

7.150256-02 

1.039406 

00 

1.055706 

00 

1.002096 

00 

1.015686 

00 

00 

1.19280F 

10“ 

“4.9  00336-  02 

17039356 

00 

1.055706 

00 

1.001946 

00 

1.01573r 

“DO 

4.937756 

00 

1. 192806 

00 

6.600376-02 

1.03  9346 

00 

1.055756 

00 

1.001776 

00 

1.015786 

00 

lU096t 

00 

1.19280E 

00 

6.  3003  fb— 02 

17039346 

■00“ 

7.7)55806 

■00“ 

- 1.00160E 

DO — 

— 1701 5846 

DO" 

5.280166 

00 

1.192806 

00 

6. 099976-02 

1.039306 

00 

1.055806 

00 

1.001506 

00 

1.0158 76 

00 

13 

0(7 

1.192796 

UU 

57899566- 02 

1.03926  EDO 

1.055806 

00 

1.001396 

00 

1.015916 

OO 

5.622576 

00 

1.192806 

00 

5.699786-02 

1.039286 

00 

1.055856 

00 

1.001306 

00 

1.015946 

00 

5.793786 

00 

1.192806 

00 

5.499996-02 

1.03929  6  DO^ 

I. 05590  E 

00 

1.001206 

ffO  ■ 

1.015986 

00 

5.964996 

00 

1. 192816 

oc 

5.375386-02 

1.039286 

00 

1.055906 

00 

1.001156 

00 

1.016006 

00 

6. 136196 

00 

1.192816 

00 

5.250776-02 

I«039826E 

■00“ 

- 1.055906 

1)0“ 

- 1.001106 

■00 — 

1.016016 

Dir 

6.307406 

00 

1.192826 

00 

5.125536-02 

1.039296 

00 

1.055956 

00 

1.00105E 

00 

1.016036 

00 

6*47  06iE 

00 

1.192826 

00 

!>•  UUU3U&-~iU 

1703932 6 

00 

1.056006 

00^ 

“17001006 

OD 

”1.014056 

00 

6.649816 

00 

1.192876 

00 

4.900316-02 

1.039306 

00 

1.056006 

00 

1.001006 

00 

1.016066 

00 

6.821U26  UO 

1.192926 

“00 

4.  800316-Q2~ 

77039296 

00 

“1.056006 

00 

1.001006 

00 

1.016086 

00 

6.992236 

00 

1.192996 

00 

4.650876-02 

1.039326 

00 

1.056056 

00 

1.001006 

00 

1.016106 

00 

7.  163436 

00 

1. 193066 

00 

“4.'5‘0140E-02 

1.039356 

w 

- 1.056106 

DOT- 

l.OOlOOE 

■00 — 

— 1.016126 

7.334646 

00 

1.193126 

00 

4.351336-02 

1.039336 

00 

1.056106 

00 

1.001006 

00 

1.016146 

00 

7.505856 

00 

T719319E“ 

00 

4.201226-02 

1.039316 

00 

1.056T06^0 

1.001006 

00 

1.016146 

00 

*M  values  are  slightly  less  than  one,  from  the  hemisphere  solution, 
but  are  adjusted  by  slightly  increasing  U  to  fit  the  requirement 
of  MOC  calculation. 
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All  M  values  are  slightly  less  than  one,  from  the  hemisphere  solution, 
but  are  adjusted  by  slightly  increasing  U  to  fit  the  requirement  of 
MOC  calculation. 
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00 

1.001006 

00 

1.0  1  4286 

00 

1.S7238E 

01 

l93dlE 

00 

27101116-02 

T.  040336  00 

17057276 

00  ■ 

1.001006 

00 

“^.0 14286 

00 

1.58950E 

01 

1. 19381E 

00 

2.10111E-02 

1.040206  00 

1.057136 

00 

1.001006 

00 

1.014286 

00 

1.60662E 

01 

I.19381E 

00 

7Z710111E-02 

17 040076  00 

17O5T00EU0 

r.ooiooFoo 

“T.OI428FOO 

1.62374E 

01 

1.19378E 

00 

2.  lOlOSE-02 

1.040136  00 

1.057006 

00 

1.001006 

00 

1.014226 

00 

1 •64U66t 

01 

1. 19375E 

OC 

2.  10099E-02 

1.040196  00 

"00“ 

- 1. 00099 E 

TTO - 

1.014176 

“DO 

1.65798E 

01 

1.19375E 

00 

2.050866-02 

1.040206  00 

1.057006 

00 

1.001006 

00 

1.014146 

00 

T167510E 

01 

r.l9375E 

00  27000736-02“ 

I704020E  00 

17057006  ^0 

““r.l701006 

00  ' 

“37014156 

00 

1.69222E 

01 

1.19382E 

00 

2.0002  56-02 

1.040046  00 

1.056756 

00 

1.001106 

00 

1.014046 

00 

1. 70934E 

01 

1. “193906 

00 

T.9997aE-02 

1.039886  00 

l«05o50c  OO 

1.001206 

00 

1.0159 86 

00 

1.72646E 

01 

1.19410E 

00 

1.949976-02 

1.039716  00 

1.056256 

00 

1.001406 

00 

1.015916 

00 

1. 74358E 

Ul 

1. 19431E 

00 

~  1.900146-02 

1^03  9546  00 

1.054016 

"00“ 

1.001406 

"DO - 

1.015846 

“00 

1.76070t 

01 

1.19450E 

00 

1.899866-02 

1.039366  00 

1.055756 

00 

1.001806 

00 

1.015776 

00 

!• JIlUZE 

01 

1.19 4 TOE  00 

T78 995  76-02 

17039176  00 

1.055506 

“00“ 

r.00l996 

“DO 

1701571F 

“DO 

1.79494E 

01 

1.19495E 

00 

1.8  748  86  -02 

1.03  902  6  00 

1.055256 

00 

1.002256 

00 

1.015426 

00 

1.81206E  01 

1.19521E 

00 

17850186-02 

r.  038876  00 

“170550  IF 

“00“ 

1.00250F770 - 

3.0I554FT70 

1.82918E 

01 

1. 19540E 

00 

1.849896-02 

1.038686  00 

1.054756 

00 

1.002706 

00 

1.015476 

00 

l»84630t 

01 

1.19560E 

00 

1.849596-02 

1.038506  OO 

— 1.054506 

170“ 

1.0O2B96 

"00 — 

1.015406 

"DO 

1.86342E 

01 

1. 1958SE 

00 

1.824876-02 

1.038346  00 

1.054256 

00 

1.003146 

00 

1.015326 

00 

I«88054E 

01 

“17196106 

00 

“17890136-02 

1.03818E  00 

1.054006 

OF 

17003396 

OO 

1.015236 

00 

1.89 766E 

01 

1.19635E 

00 

1. 799926-02 

1.03  8  03  6  00 

1.053756 

00 

1.003456 

00 

1.015156 

00 

1. 914788 

01 

1 . 1966 IE 

00 

I.  79  9  7  06  -  02 

00 

l«0535Tt  DO 

1.003906  OO 

'  1.015046 

00 

1.93190E 

01 

1.19680E 

00 

1.  77486E-02 

1.03  7696  00 

1.053256 

00 

1.004096 

00 

1.015006 

00 

1.94902E 

01 

1.19700E 

00 

1.  750016-02 

1. 03751 E  00 

— 1.05300E 

170“ 

170042SF 

"00 - 

1.014936 

“DO 

1.96614E 

01 

1.19720E 

00 

1.750316-02 

1.037346  00 

1.052756 

00 

1.004496 

00 

1.014846 

00 

uu 

im 

1.03  7176  00 

1.05251FO0  1.004706  W 

1.014796 

OO 

2.00038E 

01 

1.19765E 

00 

1.725226-02 

1.037006  00 

1.052256 

00 

1.004946 

00 

1.014716 

00 

2*  Oi  7!>0E 

01 

1.19790E 

00 

“17699826-02 

67036836  00 

1.052006 

0O“' 

17005196 

00 

“17014426 

00 

2.03463E 

01 

1.19795E 

00 

1.674736-02 

1.036856  00 

1.052006 

00 

1.005246 

00 

1.014416 

00 

2.05175E 

01 

1.19800E 

00 

1.  649656—  02 

1.034876  00 

1.052OO6 

170“ 

- 1.005286 

770 - 

1.014596 

“DO 

2.06aa7E  01 
2.(JB599E  01 


1.19007E  OC 
“1719 8r4E 


1.62458E-02 
■02“ 


1. 03691 E  00 
“r7D369«^0a 


1.05201E  00 
T705203E^011 


i.00535E  00 
I700541t^  00 


1.01457E  00 
T7Xa455E  00 


♦M  values  are  slightly  less  than  one,  from  the  hemisphere  solution, 
but  are  adjusted  by  slightly  increasing  U  to  fit  the  requirement 
of  MOC  calculation. 


AEDC-TR-75-114 


M  =  1.1 

oe 


R 

u 

V 

RHO 

p 

M 

T 

1.00000^ 

00 

1.785766 

CO 

0.0 

5.470906- 

-01 

4.300006-01 

1.702376 

00 

7.85S77E- 

-01 

oa — 

TT689286 

00 

6.784496-02 

6.331186- 

-01 

5.275056- 

-01 

1.565376 

00 

8.331856 

-01 

00 

1.597636 

90 

1.283306-01 

7.147066-01 

6.250106-01 

1.448546 

00 

8.744996 

-01 

1751 7e^P 

00  ■ 

1.509056 

00 

1.818206-01 

7.927336 

-01 

7.225156-01 

1.345586 

00 

9.114236-01 

1. 423796 

00 

1.422456 

00 

2.285166-01 

8.678086-01 

8.200206-01 

1.252586 

00 

9.449326 

-01 

l.r2974F 

00  ■ 

1.395636 

00 

2.123176-01 

8.953776 

-01 

8.566456- 

-01 

1.219776 

OC 

9.567436-01 

1.63568E 

00 

1.368946 

00 

1.965936-01 

9.226136- 

-01 

8.932706- 

-01 

1.187886 

00 

9.681956 

-01 

l.7V«3r 

"00 - 

1.342346 

"00  "" 

1.813376-01 

9.495376 

-01 

9.298956- 

-01 

1.156826 

00 

9.793146 

-01 

1.847586 

00 

1.315806 

00 

1.665406-01 

9.761616- 

-01 

9.665206-01 

1.126506 

00 

9.901246 

-01 

1.953535 

00 

1.298056 

00 

1.644076-01 

9.930306- 

-01 

9.898996-01 

1.107566 

00 

9.968476-01 

2.05947F 

90 

1.280366 

00 

1.622796-01 

1.009796 

00 

1.013286 

OC 

1.088886 

00 

1.003466 

00 

7.16942'r 

00 

1.26272" 

00 

1.601546-01 

1.026436 

00 

1.036666 

00 

1.070436 

00 

1.009966 

oc 

2.271376 

00 

1.245136 

00 

1.580326-01 

1.042986 

00 

1.060046 

00 

1.052206 

00 

1.016366 

00 

2.37732P 

"00 - 

“^7  2434Br  60  ~  " 

I. 493016-01 

1.045686 

00 

1.063796 

00 

1.049436 

00 

1.017326 

00 

2.483266 

00 

1.241776 

00 

1.405876-01 

1.048386 

00 

1.067556 

00 

1.046666 

00 

1.018296 

oc 

275S92Tf 

66 

1.240016 

oc 

1.318906-01 

1.051076 

00 

1.071316 

00 

1.043916 

00 

1.019266 

00 

2.695166 

09 

1.238186 

00 

1.232126-01 

1.053766 

00 

1.075076 

OC 

1.041156 

00 

1.020226 

00 

85 11 16 

00 

1.232236 

00 

1.210616-01 

1.059606 

00 

1.083326 

00 

1.034926 

00 

1.022396 

00 

2.90705P 

00 

1.226286 

00 

1.189246-01 

1.065426 

00 

1.091576 

00 

1.028706 

00 

1.024556 

00 

3.61350P 

"00 - 

“17225336 

00 

1. 168046-01 

1.071226 

00 

1.099836 

00 

1.022526 

00 

1.026706 

00 

3.11895'= 

00 

1.214386 

00 

1.146986-01 

1.077026 

00 

1.108086 

00 

1.016356 

00 

1.028846 

00 

3.224996 

00 

1.214986 

00 

1.099016-01 

1.077076 

00 

1.108086 

00 

1.016516 

00 

1.028796 

00 

3.339846 

00 

1.215576 

9C 

1.050986-01 

1.077136 

00 

1.108086 

00 

1.016676 

00 

1.028736 

00 

5^436796 

06 

1.216136 

00 

1.002886-01 

1.077196 

00 

1.108086 

00 

1.016836 

00 

1.028686 

00 

3. 542 74c 

00 

1.216676 

00 

9.547236-02 

1.077256 

00 

1.108086 

00 

1.016986 

00 

1.02 8626 

00 

3.64S69C 

00  ~ 

■  1.216486 

00 

9.408896-02 

1.077716 

00 

1.108676 

00 

1.016696 

00 

1.028726 

00 

3.754636 

00 

1.216296 

00 

9.270566-02 

1.078186 

00 

1.109266 

00 

1.016396 

00 

1.028836 

00 

3.960586 

90 

1.216096 

09 

9.132256-02 

1.078656 

00 

1.109866 

00 

1.016086 

00 

1.028936 

00 

3.966536 

00 

1.215896 

00 

8.993966-02 

1.079116 

00 

1.110456 

00 

1.015786 

00 

1.029046 

00 

47672486 

00 

1.214436 

00 

8.691286-02 

1.080846 

00 

1.112866 

00 

1.014096 

00 

1.029626 

00 

4. 17842E 

00 

1.212956 

00 

8.3893 76-02 

1.082566 

00 

1.115276 

00 

1.012406 

00 

1.030216 

00 

4.28437F  00  ' 

1.211476 

00 

8.097986-02 

1.084296 

00 

1.117686 

00 

1.010716 

00 

1.030796 

oc 

4.39032" 

00 

1.209976 

00 

7.787386-02 

1.086016 

00 

1.120096 

00 

1.009026 

00 

1.031386 

00 

47496276 

00 

1.210196 

00 

7.678036-02 

1.086066 

00 

1.120096 

00 

1.0C916E 

00 

1.031336 

00 

4.602216 

00 

1.210406 

00 

7.568606-92 

1.086116 

00 

1.120096 

00 

1.009306 

00 

1.031286 

00 

4.709196 

00 

1.210606 

oc 

7.459146:>02 

1.086176 

00 

1.120096 

00 

1.009456 

oc 

1.031236 

00 

4.81411" 

00 

1.210816 

00 

7.349616-02 

1.086226 

00 

1.120096 

00 

1.009586 

00 

1.031186 

00 

4.920066 

"00 — 

“17211066 

00 

7.168486-02 

1.086276 

00 

1.120096 

00 

1.009726 

00 

1.031146 

00 

5.02600" 

00 

1.211306 

00 

6.987366-02 

1.086326 

00 

1.120096 

00 

1.009866 

00 

1.031C96 

00 

♦Mach  numbers  are  slightly  less  than  one,  from  the  hemisphere  solution, 
but  are  adjusted  by  slightly  increasing  U  to  fit  the  requirement  of 
MOC  calculation. 


AEDC-TR-75-114 


M  =  1.1  (Continued) 


R 

U 

V 

RHO 

p 

M 

T 

3713 195P 

00 

1.21153? 

00 

6.806036-02 

1.08636? 

00 

1.12009E 

00 

1.00999? 

00 

1.03104? 

OC 

5.23790E 

00 

1.211776 

00 

6.624736-02 

1.08641E 

00 

1.12009E 

00 

1.01012? 

00 

1.03100? 

00 

3;343e5P 

00  - 

~  1.21195? 

00 

6.536286-02 

1.08646? 

00 

1.12009? 

00 

1.01026? 

00 

1.03095? 

00 

5.44979E 

00 

1.21212? 

00 

6.44788E-02 

1.086516 

00 

1.12009E 

00 

1.01038? 

00 

1.03091? 

00 

3753373r 

00 

“  T7212306 

00 

6.359356-02 

1.98655E 

00 

1.12009E 

00 

1.01051? 

00 

1.03086? 

00 

5.66169E 

00 

1.21247? 

00 

6.270906-02 

1.08660E 

00 

1.12009E 

00 

1.01064? 

00 

1.03C82E 

OC 

5775764F' 00 

1.21316F 

00 

6.153136-02 

1.08619? 

00 

1.11943? 

00 

1.011266 

00 

1.03C60? 

00 

5.87358'; 

00 

1.21384? 

00 

6.035216-02 

1.08578? 

00 

1.11877E 

00 

1.01189? 

00 

1.03039? 

00 

3737353? 

DO 

1.21452? 

00 

5.917146-02 

1.08536? 

00 

l.liailE 

00 

1.01251? 

00 

1.03017? 

00 

6. 0854 BE 

00 

1.21520? 

00 

5.79892E-02 

1.08495? 

00 

1.11745E 

OC 

1.013136 

00 

1.02995? 

00 

37T3I435 

■W 

1.21*485? 

00 

5;72377E-02 

1.08545E 

00 

1.11811? 

00 

1.01275? 

00 

1.03CC9? 

00 

6.29737E 

00 

1.21450? 

00 

5.648656-02 

1.08595E 

00 

1.11877? 

00 

1.01237? 

00 

1.03022? 

00 

374333 

DO  “ 

1.21615E 

00 

5.573566-02 

1.08645E 

00 

1.11943E 

00 

1.011986 

00 

1.03C35? 

00 

6.50927E 

00 

1.21380? 

00 

5.49850E-02 

1.08695E 

00 

1.12009? 

00 

1.01159? 

00 

1.03C49? 

00 

6.61 522E 

00 

1.21395? 

00 

5.417246-02 

1.08699? 

00 

r.r2do9? 

00 

1.01171? 

00 

1.03045? 

00 

6.72U6E._30 

1.214106 

00 

5.33595E-02 

It08703? 

00 

'  1,12009E 

OQ 

_  UQil82E 

00 

1.03041? 

00 

6.82711C 

00 

1.21624? 

00 

5.254626-02 

1.08707E 

00 

1.12009E 

00 

1.01193? 

00 

1.03037? 

00 

6.93306E 

00 

1.216386 

00 

5. 17326? -02 

1.08711? 

00 

1.12009? 

00 

1.012036 

00 

1.03034? 

00 

7703 96 IE" 

00 

' 1.21451? 

00 

5. 11249? -02 

1.08715? 

00 

1.12009? 

00 

1.01214? 

00 

1.03030? 

00 

7.14495= 

00 

1.21664? 

00 

5.05170E-02 

1.087 19E 

00 

1.12009? 

00 

1.01225? 

00 

1.03026? 

00 

T.  250906 

00 

1.21477? 

00 

4.990886-02 

1.0e722E 

00 

1.12009E 

00 

1.01235? 

00 

1.03023? 

00 

7.35685E 

00 

1.21490? 

00 

4.930056-02 

1.08726E 

00 

1.12009? 

00 

1.01245? 

00 

1.03019? 

00 

7.462803  00 

— I.215D2? 

00  — 

4.871863-02 

■  1.0DT30? 

00 

1.12009? 

00 

1.01255? 

00 

1.03016? 

00 

7.56874? 

00 

1.21514? 

00 

4.813786-02 

1.08733? 

00 

1.12009? 

00 

1.01265? 

00 

1.03012? 

00 

7.67663S 

DO' 

1.21526? 

00 

4.75555E-02 

1. 087376 

00 

1.12009? 

00 

1.01274? 

00 

1.03009? 

00 

7.780646 

00 

1.21538? 

00 

4.697446-02 

1.08740E 

00 

1.12009E 

00 

1.01284? 

00 

1.030066 

00 

7.88659r 

00 

1.21599? 

00 

4.64934E-02 

1.08698? 

00 

1.11943? 

00 

1.013446 

00 

1.02985? 

00 

7.99253? 

00 

1.21660? 

OC 

4.60119E-02 

1.08656? 

00 

1.11877? 

00 

1.01403? 

00 

1.02965? 

00 

8. 0984 ST 

“00“ 

1.21721? 

00 

4.552976-02 

1.08613? 

00 

1.11811? 

00 

1.01462? 

00 

1.02944? 

00 

8.20443? 

00 

1.21781? 

00 

4.50470E-02 

1.08571E 

00 

1.11745E 

00 

1.01521? 

00 

1.02924? 

00 

«73n}38f 

scr- 

1.21799? 

00 

4.46173E-02 

1.08567E 

00 

1.11736E 

00 

1.01537? 

00 

1.02918? 

OC 

8.41632? 

00 

1.21815? 

00 

4.418866-02 

1.08565? 

00 

1.11727E 

00 

1.01552? 

00 

1.02913? 

00 

f732227?~00 

1.21832? 

00 

4.375866-02 

1.08561? 

00 

1.11718E 

00 

1.01568? 

00 

1.02908? 

OC 

8.62822? 

00 

1.218496 

00 

4.332956-02 

1.08558? 

00 

1.11709E 

00 

1.01583? 

00 

1.029026 

00 

8.73417E 

“OXT" 

T."210333“Dir 

4.29150E-02 

*1.08561? 

00 

■“  "  1.11709? 

00 

1.01591? 

00” 

— 1.02900? 

00 

8.84011? 

00 

1.21868? 

00 

4.25015E-02 

1.08564E 

00 

1.11709? 

00 

1.01599? 

00 

1.02897? 

00 

37343063 

00  “ 

1.2l877r  OC 

4.208666-02 

1.08567? 

00 

1.11709? 

00 

1.01607? 

00 

1.02894? 

00 

9.05201? 

00 

1.21886? 

00 

4.167296-02 

1.08569? 

00 

1.11709? 

OC 

1.01614? 

00 

1.02891? 

00 

37T5796T  06 

1.21933? 

00 

4.134986-02 

1.08538? 

00 

1.1 1659? 

00 

1.01660? 

00 

1.02876? 

CO 

♦Mach 

numbers  are 

slightly  less 

than  one 

9 

from  the 

hemisphere 

solution, 

but  are  adjusted  by  slightly  increasing  U  to  fit  the  requirement  of 
MOC  calculation. 


AEDC-TR-75-114 


M  =1.1  (Concluded) 


R 

U 

V 

RHO 

p 

M 

T 

9.26390P 

00 

1.21979?  00 

4.10276E-02 

1.08506? 

00 

1.11609? 

00 

1.01705? 

00 

1.02860? 

00 

9.36985? 

00 

1.22025?  00 

4. 0703 9? -02 

1.08473? 

00 

1.11559? 

00 

1.01750? 

00 

1.02844? 

00 

9.47580C 

00 

1.22071?  00 

4.03810E-02 

1.08441? 

00 

1.11509? 

00 

1.01795? 

00 

1.02829? 

00 

9.58175? 

HO 
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APPENDIX  B 

ADDITIONAL  DATA  FOR  =  1.05,  1.1,  1.2,  AND  1.5 

Because  AEDC  data  were  consistently  lower  in  shock  standoff  distance  than 
predicted  by  theory  in  the  low  supersonic  flow  regime,  additional  shadowgraph  and  pressure 
data  were  obtained  with  the  same  testing  setup  in  the  Aerodynamic  Wind  Tunnel  ( I T)  for 
=  1 .05,  1.1 , 1.2,  and  1 .5  after  the  completion  of  the  technical  report.  Therefore,  the  new 
data  are  included  herein  as  an  appendix. 

In  the  additional  experiments  special  attention  was  given  to  the  optical  lineup  in  taking 
shadowgraphs,  and  the  area  covered  by  the  pictures  was  increased  to  1 1  by  14  in.  Figure  B-1 
shows  the  shadowgraphs. 

A  comparison  of  the  present  data  with  the  theoretical  results  reported  in  the  main  text 
is  given  in  Fig.  B-2  for  the  shock  standoff  distance  and  in  Fig.  B-3  for  the.shock  position. 
Two  sets  of  .data,  one  from  each  side  of  the  model  at  different  runs,  were  included  for  M  = 
1 .2  only.  It  is  seen  that  a  better  agreement  between  theory  and  experiment  is  obtained.  No 
significant  change  in  the  surface  pressure  between  the  more  recent  data  and  the  data  shown 
in  Fig.  2  for  M^=  1.05,  1.1,  and  1.2  can  be  found.  The  more  recent  data,  therefore,  provide 
further  evidence  to  confirm  the  conclusion  presented  in  the  main  text. 
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Figure  B-3.  Comparison  of  theoretical  shock  position 
with  experiment,  7  =  1.4. 
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NOMENCLATURE 

a  Local  speed  of  sound 

Cp  Surface  pressure  coefficient 

Cp  Specific  heat  of  the  gas  at  constant  pressure 

Cy  Specific  heat  of  the  gas  at  constant  volume 

K  Curvature  of  the  body  in  the  meridian  plane 

M  Mach  number 

n,s  Body-oriented  coordinate  system 

P  Pressure 

q  Total  velocity 

R  Radius  of  the  cylinder 

r,Y  Radial  distance  from  the  body  axis 

5  Entropy  of  the  gas 

u,v  Velocity  components  in  the  n  and  s  directions,  respectively 

Z  Axial  distance  from  the  nosetip 

7  Specific  heat  ratio,  Cp/cy 

A  Bow  shock  standoff  distance 

* 

Tj,^  Coordinates  in  the  computation  plane 

6  Inclination  angle  of  the  body 

K  1  +  Kn 

H  Mach  number 

p  Density  of  the  gas 

f}/  Velocity  potential 

(  Subscript  for  free-stream  conditions 
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